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Abstract. Given a matrix A G F(t)[D; S] nxn over the ring of differential 
polynomials, we show how to compute the Hermite form H of A and a 
unimodular matrix U such that UA — H . The algorithm requires a 
polynomial number of operations in F in terms of n, deg^,^, deg t A 
When F = Q it require time polynomial in the bit-length of the rational 
coefficients as well. 



1 Introduction 

Canonical forms of matrices over principal ideal domains (such as Z or F[x], 
for a field F) have proven invaluable for both mathematical and computational 
purposes. One of the successes of computer algebra over the past three decades 
has been the development of fast algorithms for computing t hese canonical fo rms. 
These include triangular forms su ch as the Her mite form (jHermitel . 1863 ). low 
degre e forms like t he Popov form (jPopovl . Il972n . as well as the diagonal Smith 
form (|Smithl . ll86llh 

Canonical forms of matrices over non-commutative domains, especially rings 
of differential and difference operators, are also extreme l y usef ul. T hese have 
been e xamined at least since Dickson (1923), Wedderburnl (1932), and ljacobsonl 
( 1943f ). A typical domain under consideration is that of differential polynomials. 
For our purposes these are polynomials over a function field F(t) (where F is a 
field of characteristic zero, typically an extension of Q, or some representation 
of C). A differential indeterminate D is adjoined to form the ring of differential 
polynomials F(t)[T>;d], which consists of the polynomials in F(t)[2?] under the 
usual addition and a non-commutative multiplication defined such that T>a = 
aV + S(a), for any a G F(t). Here S : F(t) — > F(t) is a pseudo-derivative, a function 
such that for all a, b G F(t) we have 



S(a + b) = 6(a) + 6{b) and S(ab) = aS(b) + 5(a)b. 

The most common derivation in F(t) takes 5(a) = a' for any a G F(f), the usual 
derivative of a, though other derivations (say S(t) = t) are certainly of interest. 

A primary motivation in the definition of F(t)[.D;5] is that there is a nat- 
ural action on the space of infinitely differentiable functions in t, namely the 
differential polynomial 



V 7 ' 



aiP + ao G F(t)[V;5] 
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acts as the linear differential operator 



a m (t) 



d m y(t) 
dt m 



l y(t) 



dt r ' 



ai(t) 



dy(t) 
dt 



a (t)y(t) 



on a differentiable function y(t). Solving and analyzing systems of such opera- 
tors involves working with matrices over F(£)[2?;<5], and invariants such as the 
differential analogues of the Smith, Popov and Hermite forms provide important 
structural information. 

In commutative domains such as Z and F[x], it has been more common to 
compute the triangular Hermite and diagonal Smith form (as well as the lower 
degree Popov form, especially as an intermediate computation). Indeed, these 
forms are more canonical in the sense of being canonical in their class under mul- 
tiplication by unimodular matrices. Polynomial-t ime algor i thms for the Smith 
and Hermite for ms over F\x] were de v eloped b y Kannanl (1985f). with impor- 
tant advances bv lKaltofen et all (ll987MVillardl ijlQQ.^ .I Mulders and Storjohannl 
and many others. One of the key features of this recent work in com- 
puting normal forms has been a careful analysis of the complexity in terms of 
matrix size, entry degree, and coefficient swell. Clearly identifying and analyz- 
ing the cost in terms of all these parameters has led to a dramatic drop in both 
theoretical and practical complexity. 

Computing the classical Smith and Hermite forms of matrices over differ- 
ential (and more general Ore) domains has received less attention though nor- 
mal forms of differential polynomial matrices have applications in so lving dif- 
ferential systems and control theory. EE ramov and Bronsteinl (|200lh analyzes 
the number of red u ction steps necessary to compute a row-reduced form, while 
Beckermann et al. ( 20061 ) analyze the complexity of row reduction in terms of 



matrix size, degree and the sizes of the coefficients of some shifts of the input 
matrix. iBeckermann et al I (120061) demonstrates tight bounds on the degree and 
coefficient siz es of the output, which we will employ here. For the Popov form, 
Chend (|2003h gives an algorithm for matrices of shift polynomials. Cheng's ap- 



proach involves order bases computation in order to elimina t e lowe r order terms 
of Ore polynomial matrices. A main contribution of [Chen g (2003) is to give an 
algorithm computing the row rank and a row-reduced basis of the left nullspace 
of a matrix of Ore polynomials in a fraction-free way. This idea is extended in 
iDavies et al.l (120081) to compute Popov form of general Ore polynomial matrices. 



( 20081 ). they reduce the problem of computing Popov form to a 



nullspace computation. However, though Popov form is useful for rewriting high 
order terms with respect to low order terms, we want a different normal form 
more suited to solving system of linear diophantinc equations. Since the Hermite 
form is upper triangular it meets this goal nicely, not to ment ion the fa c t that 
it is a "classical" canonical form. In a slightly different vein, iMiddekd (|2008) 
has recently given an algorithm for the Smith (diagonal) form of a matrix of 
differential polynomials, which requires time polynomial in the matrix size and 
degree (but the coefficient size is not analyzed). 

In this paper, we first discuss some basic operations with polynomials in 
F(t)[D; S], which are typically written with respect to the differential variable V 
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as 



/ = fa + faV + faV 2 + • • • + faV\ 



(1.1) 



where fa, . . . , fa S F(t), with fa ^ 0. We write d = deg-p/ to mean the degree in 
the differential variable, and generally refer to this as the degree of /. Since this is 
a non-commutative ring, it is important to set a standard notation in which the 
coefficients fa, . . . , fa £ F(t) are written to the left of the differential variable T>. 
For u, v £ F[t] relatively prime, we can define deg t (u/v) — max{deg t w, deg t w}. 
This is extended to f £ F(t)[D;5] as in (jl.l|) by letting deg t / = maxi{deg t /i}. 
We think of deg t as measuring coefficient size or height. Indeed, with a little 
extra work the bounds and algorithms in this paper are effective over Q(t) as 
well, where we also include the bit-length of rational coefficients, as well as the 
degree in i, in our analyses. 

A matrix U £ F(t)[D;6] n " is said to be unimodular if there exists & V £ 
F(t)[T>; 5} n ™ such that UV = I, the n X n identity matrix. Note that we do not 
employ the typical determinantal definition of a unimodular matrix, as there 
is no easy notion of determinant for matrices over F(t)[D]$\ (indeed, working 
around this deficiency suffuses much of our work). 

A matrix H £ F(t)[D;S] nxn is said to be in Hermite form if H is upper 
triangular, if every diagonal entry is monic, and every off-diagonal entry has 
degree less than the diagonal entry below it. As an example, the matrix 



Note that the Hermite form may have denominators in t. Also, while this example 
does not demonstrate it, it is common that the degrees in the Hermite form, in 
both t an T>, are substantially larger than in the input. 

In this paper we will only concern ourselves with matrices in F(t)[T>;S] n n 
of full row rank, that is, matrices whose rows are F(t)[D; <5]-linear independent. 
For any matrix A £ F(t)[D;S\ ™, we show there exists a unimodular matrix 
U such that UA = H is in Hermite form. This form is canonical in the sense 
that if two matrices A, B £ F(t)[D; S] nxn are such that A = PB for unimodular 
P £ F(t)[T>; S] n n then the Hermite form of A equals the Hermite form of B. 

The main contribution of this paper is an algorithm that, given a matrix 
A £ F(t)[D;S] n n (of full row rank), computes H and U such that UA = H, 
which requires a polynomial number of F-operations in n, deg-pA, and Aeg t A. It 
will also require time polynomial in the coefficient bit-length when F = Q. 

The remainder of the paper is organized as follows. In Section 2 we summarize 
some basic properties of differential polynomial rings and present and analyze 




has Hermite form 



2 + t + V 1 + 2* 
2 + t + V 





4 



algorithms for some necessary basic operations. In Section 3 we introduce a new 
approach to compute appropriate degree bounds on the coefficients of H and 
U. In Section 4 we present our algorithm for computing the Hermite form of a 
matrix of differential polynomials and analyze it completely. 



2 Basic structure and operations in F[t][X>;(5] 

In this section we discuss some of the basic structure of the ring F(t)[2?;5] and 
present and analyze simple algorithms to do some computations that will be 
necessary in the next section. 

Some well-known propert ies of F(t)[D;5] are worth recalling; see 



Bronstein and Petkovsekl (|1994l ) for an algorithmic presentation of this theory. 



Given f,g£ F(t)[D; 8] , there is a degree function (in T>) which satisfies the usual 
properties: deg v (fg) = deg v f + deg v g and deg v {f + g) < max{deg£,/, deg v g}. 
F(t)[T>; 5} is also a left and right principal ideal ring, which implies the existence 
of a right (and left) division with remainder algorithm such that there exists 
unique q,r G F(t)[V;5] such that / = qg + r where deg-p(r) < deg v (g). This 
allows for a right (and left) euclidean-likc algorithm which shows the existence of 
a greatest common right divisor, h — gcrd(/, g), a polynomial of minimal degree 
(in T>) such that / = uh and g = vh for u, v G F(t)[D; 8]. The GCRD is unique 
up to a left multiple in F(i)\{0}, and there exist co-factors a,b G F(t)[X>;<5] 
such that af + bg = gcrd(/, g). There also exists a least common left multiple 
lclm(/, g). Analogously there exists a greatest common left divisor, gcld(/, g), 
and least common right multiple, lcrm(/, g), both of which are unique up to a 
right multiple in F(t). 

Efficient algorithm s for computing products of polynomials are developed in 
van der Hoevenl (|2002l ) and iBostan et al. (l2008h. while fas t algo rithms tocom- 



pute the LCLM and GCRD, are developed in lLi and Nemed (|l997l ) andO (|l998l ). 
In this paper we will only need to compute very specific products of the form 
T> k f for some k G N. We will work with differential polynomials in F[i][£>; S], as 
opposed to F(t)[T>;5), and manage denominators separately. If / G F[f][2?;<$] is 
written as in (jl.ip . then / , . . . , fd G F[t], and 

0<i<d 0<i<d 

where /■ G F[t] is the usual derivative of fi G F[t}. Assume deg t / < e. It is easily 
seen that degp(2?/) = d + 1, and deg t (Vf) < e. The cost of computing Vf is 
O(de) operations in F. Computing T> k f , for 1 < k < m then requires 0{dem) 
operations in F. 

If F = Q we must account for the bit-length of the coefficients as well. Assum- 
ing our polynomials are in Z[t][2?;i5] (which will be sufficient), and are written 

as above, we have /; = J2o<j< e fifi for fa e Z - We write 11/11°° = max \fv I to 
capture the coefficient size of /. It easily follows that HP/Hoo < (e + l)||/||oo, 
andso ||P m /l|oo < (e + ini/lloo- 
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Lemma 2.1. 

(i) Let f G F[i][X>;<5] have deg v f = d, dcg f / = e, and let m G N. Then we can 

compute D k f , for 1 < k < m, with O(dem) operations in F. 
(ii) Let f G Z[t][V;5}. Then ||X>"7lloo < (e + l) m ■ \\f\\oo, and we can compute 
T) 1 f , for 1 < i < m, with 0(dem ■ (mloge + log ||/||oo) 2 ) bit operations. 

We make no claim that the above methods are the most efficient, and faster 
polynomial and matrix arithmetic will certainly improve the cost. However, the 
above analysis will be sufficient, and these costs will be dominated by others in 
the algorithms of later sections. 



3 Existence and degree bounds on the Hermite form 

In this section we prove the existence and uniqueness of the Hermite form over 
F(t)[D; (5], and prove some important properties about unimodular matrices and 
equivalence over this ring. The principal technical difficulty is that there is no 
natural determinant function with the properties found in commutative linear 
algebra. The determinant is one of the main tools used in the analysis of essen- 
tially all fast algorithms for computing the Hermite form H and transformation 
matrix U , and sp ecifi cally two relevant tec hniques in established methods by 



matr ix U , and sp ecificall y two relevant tec hniques m established methods by 
IStorkffiarnil (|l994h and lKaltofen et all (Il987t). On e appr oach might be to employ 



the non-commutative determinant of lDieudonni (1943), but this adds consider 



able complication. Instead, we find degree bounds via established bounds on the 
row-reduced form. 



Definition 3.1 (Unimodular matrix). Let U G F(t)[T>;5] n n and suppose 
there exists a V G F(t)[D; 6] n n such that UV — L n , where L n is the identity 
matrix over F(t)[D; 5} nxn . Then U is called a unimodular matrix over F(t)[D; 5]. 

This definition is in fact symmetric, in that V is also unimodular, as shown in 
the following lemma (the proof of which is left to the reader). 



Lemma 3.1. Let U G F(t)[D; 5} n n be unimodular such that there exists a V G 
F(t)[V;S] nxn with UV = I n . Then VU = I n as well. 



Theorem 3.1. Let a,b G F(i)[Z?;<5]. There exists a unimodular matrix 

W=(^ s V ^eF{t)[V-5] 2y2 suchthat W Q = 
where g = gcrd(a, b) and sa = —tb = lclm(a, b). 

Proof. Let u, v G F(t)[T>; 5] be the multipliers from the euclidean algorithm such 
that ua + vb = g. Since sa = —tb = lclm(a,6), we know that gcld(s,i) = 1 
(otherwise the minimality of the degree of the lclm would be violated). It follows 
that there exist c, d G F(t)[T>; 5} such that sc + td = 1. Now observe that 

(u v\ ( ag^ 1 c\ fl —uc — vd\ (\ uc + vd\ ( 1 — uc — vd\ / 1 0\ 
I s t [bg- 1 d)\0 1 J = { 1 JU 1 J = U I ■ 
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Thus 



w- 1 = 



ag 1 ag 1 {—uc — vd)-\-c 
bg^ 1 bg^ 1 (— ttc — vd) + d 



ag —a 
bg- 1 -b- 



so W is unimodular. 



Definition 3.2 (Hermite Normal Form). Let H G F(t)[V; S] nxn with full 
row rank. The matrix H is in Hermite form if H is upper triangular, if every 
diagonal entry of H is monic, and if every off- diagonal entry of H has degree 
(in T>) strictly lower than the degree of the diagonal entry below it. 



Let A G 

nxn 



Theorem 3.2 

matrix H G F(t)[V;S] 
matrix U € F(i)[£>;,5] ? 



F(t)fD; S] nxn have row rank n. Then there exists a 
with row rank n in Hermite form, and a unimodular 
such that UA = H . 



Proof. We show this induction on n. The base case, n = 1, is trivial and we 
suppose that the theorem holds for n-lxn-1 matrices. Since A has row rank 
n, we can find a permutation of the rows of A such that every principal minor of 
A has full row rank. Since this permutation is a unimodular transformation of 
A, we assume this property about A. Thus, by the induction hypothesis, there 



exists a unimodular matrix U% G F(i)[2?;J] 



(n-l)x(n-l) 



such that 











A = H = 



\0 ••■01/ 



H 



2,2 





V^n.l ^n,2 



#7i-l,n-l * 
-^n,n— 1 A n ^ n J 



€F(t)[D;S] 



where the (n — l)st principal minor of H is in Hermite form. By Theorem 13.1 
wc know that there exists a unimodular matrix 



W = 



-t, 



G F(t)[V;5] 2x2 such that W 



Hi; 

A n ,_ 



€F(t)[V;6] 



This allows us to reduce A n> i, . . . ,^4 n! „_i to zero, and does not introduce any 
non-zero entries below the diagonal. Also, all off-diagonal entries can be reduced 
using unimodular operations modulo the diagonal entry, putting the matrix into 
Hermite form. □ 

Corollary 3.1. Let A G F(t)[D;S] have full row rank. Suppose UA — H for 



unimodular U G F(t)[T>;5] 
U and H are unique. 



and Hermite form H G F(t)[T>; 6] n . Then both 



Proof. Suppose H and G are both Hermite forms of A. Thus, there exist uni- 
modular matrices U and V such that UA = H and VA = G, and G = WH 
where W = VU^ 1 is unimodular. Since G and H are upper triangular matrices, 
we know W is as well. Moreover, since G and H have monic diagonal entries, 
the diagonal entries of W equal 1. We now prove W is the identity matrix. By 
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way of contradiction, first assume that W is not the identity, so there exists an 
entry Wij which is the first nonzero off-diagonal entry on the ith row of W. 
Since i < j and since Wu = 1, Gy = + WijHjj. Because Wij ^ 0, we see 
deg v Gij > deg v Gjj, which contradicts the definition of the Hermite form. The 
uniqueness of U follows similarly. □ 

Definition 3.3 (Row Degree). A matrix T G F(t)[D; 5] nxn has row degree 
u G (N U {— oo})" if the ith row of T has degree u,-. We write rowdeg u . 

Definition 3.4 (Leading Row Coefficient Matrix). Let T G F(t)[D;S] nxn 
have rowdegu. Set N = deg v T and S = diag(V N ~ Ul , . . . , V N ~ U "). We write 

ST = LT> N + lower degree terms in T>, 

where the matrix L = LC row (T) G F(t) nxn is called the leading row coefficient 
matrix of T. 

Definition 3.5 (Row-reduced Form). A matrixT G F(t)[D; S] mxs with rank 
r is in row-reduced form if rank LC row (T) — r. 



Fact 3.1 (jBeckermann et al . (2006) Theorem 2.2). For any A G F(t)[D;S] 7 



there exists a unimodular matrix U G F(i)[P; S] mxm , with T = UA having 
r < min{r7i, s} nonzero rows, rowdegT < rowdeg A, and where the submatrix 
consisting of the r nonzero rows of T are row-reduced. Moreover, the unimodular 
multiplier satisfies the degree bound 



rowdegU < V + (\Tt\ — IT?] — min{uj}) 



e 



where ~u :— max( , rowdeg A), ~v* :— max( , rowdegT) , and ~e is the column 
vector with all entries equal to 1. 

The proof of the following is left to the reader. 

Corollary 3.2. If A £ F(t)[D; 5} nxn is a unimodular matrix then the row re- 
duced form of A is an identity matrix. 

The following theorems provide degree bounds on H and U. We first compute 
a degree bound of the inve rse of U by using the ide a of backward substitution, 
and then use the result of B eckermann et al. I (|2006h to compute degree bound 
of U. 

Theorem 3.3. Let A G F(t)[D;S] nxn be a matrix with degpAy < d and full 
row rank. Suppose UA = H for unimodular matrix U G F(t)[D;6] nxn and 
H G F(t)[X>;<5] " in Hermite form. Then there exist a unimodular matrix 
V G F(i)[D; 6} nxn such that A = VH where UV = I n and deg^Vy < d. 

Proof. We prove by induction on n. The base case is n = 1. Since H\\ = 
gcrd(An, . . . ,A n i), deg-p-ffn < d and so deg-pVn < d for 1 < i < n. Now, 
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we suppose that our claim is true for k where 1 < k < n. Then we have to show 
that deg v Vik+i < d. We need to consider two cases: 
Case 1: deg v Vi^+i > max(deg P Ki, • ■ ■ , deg v Vik)- Since 

deg v Hk+i,k+i > max(deg X) H 1)k+ x, . . . ,deg v H kik+1 ), 
degj^fe+i = degj,{V if k+\Hk+i,k+i), 

where A iik +i = VnH 1:k+1 H h Vi, k+ \H k+l , k+l . Thus, deg^K.k+i < d. 

Case 2: deg-pV^+i < max(deg- D l^i, . . . , degpVifc). Thus, by induction hypothe- 
sis, degpVi.fc+i < d. □ 

Corollary 3.3. Let A, V, and U be those in Theorem \ 3.3[ Then degpt/jj < 
(n-l)d. 

Proof. By Corollary |3.2l we know that the row reduced form of V isl n . Moreover, 
since /„ = UV, we can compute the degree bound of U by using Fact 13. II Clearly, 

V + (\~u\ — \~v\ — mhv{uj})~e < V + (\u \ — mm{uj})~e , 

j j 

where ~u := max( , rowdegV) and ~v* := max( , rowdegI n ) = . Since the 
degree of each row of V is bounded by d, (\u \ — minj{uj}) < (n ~ l)d. Then, 
by Fact 13.11 rowdegU < (n — l)d. Therefore, deg-pt/y < (n — l)d. □ 

Corollary 3.4. Let H be same as that in Theorem \3.3\ Then deg v Hij < nd. 

Proof. Since deg v Uij < (n — \)d and deg v Aij < d, deg v Hij < nd. □ 



4 Computing Hermite forms by linear systems over F(i) 

In this section we present our polynomial-time algorithm to compute the Hermite 
form of a ma trix over F(t) [D ; S] ■ W e ex hibit a variant of th e linear system method 
developed in lKaltofen et al.1 ( 1987t ) and Storjohann ( 19941 ). The approach of these 



papers is to reduce the problem of computing the Hermite of matrices with 
(usual) polynomial entries in F[z] to the problem of solving a linear system 
equations over F. Analogously, we reduce the problem of computing the Hermite 
form over F[t] [£>; 5} to solving linear systems over F(t). The point is that the field 
F(t) over which we solve is the usual, commutative, field of rational functions. 

For convenience, we assume that our matrix is over F[f][2?;5] instead of 
F(t)[T>; 5], which can easily be achieved by clearing denominators with a "scalar" 
multiple from F[i\. This is clearly a unimodular operation in the class of matrices 
over F(t)[V;S]. 

We first consider formulating the computation of the Hermite form a matrix 
over F(t)[D; S] as the solution of a "pseudo" -linear system over F(t)[T>; 5} (i.e., a 
matrix equation over the non-commutative ring F(t)[T>; 5]). 

Theorem 4.1. Let A e F[t] [D; 5} nxn have full row rank, with deg^Aij < d, 
and (di, . . . , d n ) 6 N™ be given. Consider the system of equations PA = G, for 
n x n matrices for P, G £ F(t)[D; 5] restricted as follows: 



9 



— The degree (inT>) of each entry of P is bounded by (n — l)d + maxi<i<„ di. 

— The matrix G is upper triangular, where every diagonal entry is monic and 
the degree of each off-diagonal entry is less than the degree of the diagonal 
entry below it. 

— The degree of the ith diagonal entry of G is di. 

Let H be the Hermite form of A and (hi, . . . , h n ) £ N n be the degrees of the 
diagonal entries of H . Then the following are true: 

(a) There exists at least one pair P,G as above with PA — G if and only if 
di > hi for 1 < i < n. 

(b) If di = hi for 1 < % < n then G is the Hermite form of A and P is a 
unimodular matrix. 



Proof. The proof is similar to that of iKaltofen et alJ (|1987I ). Lemma 2.1. Given 



a degree vector (d\, . . . , d n ), we view PA — G as a system of equations in the 
unknown entries of P and G. Since H is the Hermite form of A, there exist a 
unimodular matrix U such that U A = H . Thus PU~ 1 H = G and the matrix 
PXJ~ X must be upper triangular since the matrices H and G are upper triangular. 
Moreover, since the matrix PU^ 1 is in F(t)[D; S] n n , and Ga = (PU~ 1 )u • Ha 
for 1 < i < n, we know di > hi for 1 < i < n. For the other direction, we 
suppose di > hi for 1 < i < n. Let D = diag(2? dl ~' 11 , . . . ,V d "~ h "). Then since 
(DU)A = (DH), we can set P = DU and G = DH as a solution to PA = G, 
and the ith diagonal of G has degree di by construction. By Corollary 13. 3i we 
know deg-pC/ij < (n — l)d and so deg-pP^.j < (n — l)d + maxi<i<„ di. 

To prove (b), suppose di = hi for 1 < i < n and that, contrarily, G is not the 
Hermite form of A. Since PU^ 1 is an upper triangular matrix with ones on the 
diagonal, PU^ 1 is a unimodular matrix. Thus P is a unimodular matrix and, 
by Corollary 13. 11 G is the (unique) Hermite form of A, a contradiction. □ 



Lemma 4.1. Let A, P, (di,...,d n ) 7 and G be as in Theorem \4-l\ and let 
[3 := (n — l)c? + maxi<i<„ d;. Also, assume that deg t Aij < e for 1 < i,j < n. 
Then we can express the system PA = G as a linear system over F(t) as PA = G 
where 

P £ f(t) nxn( - l3+1 \ A G f^ n (P+ 1 ) xn (f ) + d + 1 ) j Q £ pu\nxn{0+d+l) _ 

Assuming the entries A are known while the entries of P and G are indeter- 
minates, the system of eguations from PA — G for the entries of P and G is 
linear over F(t) in its unknowns, and the number of equations and unknowns is 
0(n 3 d). The entries in A are in F[t] and have degree at most e. 

Proof. Since deg^P^j < /3, each entry of P has at most (/3+1) coefficients in F(t) 
and can be written as P l} = £o<fc</3 P l3 kD k . We let P e F(t) nxn ^ +1 ^ be the 
matrix formed from P with P^ replaced by the row vector (Pijo, . . . , Pijp) & F(i). 

Since deg^P < /3, when forming PA, the entries in A are multiplied by Tr for 
< £ < P, resulting in polynomials of degree in T> of degree at most ji = (3 + d. 
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Thus, we construct A as the matrix formed from A with Aij replaced by the 
(P + 1) x (a + 1) matrix whose £th row is 

(4}o^,-,Ag) such that V^A^A^+Agv + .-. + Agv^. 

Note that by Lemma HOI we can compute T> e Aij quickly. 

Finally, we construct the matrix G. Each entry of G has degree in T> of degree 
at most nd < n(f3 + d + 1). Thus, initially G is the matrix formed by G with Gjj 
replaced by 

(G ij0 , G ljtl ) where G i} = G ij0 + G vjl V + h G.^V. 

However, because of the structure of the system we can fix values of many of 
the entries of G as follows. First, since every diagonal entry of the Hermite form 
is monic, we know the corresponding entry in G is 1. Also, by Corollary 13.41 
the degree in T> of every diagonal entry of H is bounded by nd, and every off- 
diagonal has degree in T> less than that of the diagonal below it (and hence less 
than nd), and we can set all coefficients of larger powers of T> to in G. 

The resulting system PA — G, restricted as above according to Theorem l4.1[ 
has 0(n 3 d) linear equations in 0(n 3 d) unknowns. Since the coefficients in A are 
all of the form T> e Aij, and since this does not affect their degree in t, the degree 
in t of entries of A is the same as that of A, namely e. □ 

With more work, we believe the dimension of the syste m can be r e duced 
to 0(n 2 d) x 0(n 2 d) if we apply the techniques presented in Storiohannl ( 1994 ) 



Section 4.3, wherein the unkno wn coefficients of G are removed from the system. 
See also Labhalla et al. ( 19961 ). 



So far, we have shown how to convert the differential system over F(t)\D; S] 
into a linear system over F(t). Also, we note, by Theorem 14. 1[ that the correct 
degree of the ith diagonal entry in the Hermite form of A can be found by 
seeking the smallest non-negative integer k such that PA = G is consistent 
when deg-pGjj = nd for j = 1, . . . , i — 1, i + 1, . . . , n and k < deg-pG^. Using 
binary search, we can find the correct degrees of all diagonal entries by solving 
at most 0(n log (nci)) systems. We then find the correct degrees of the diagonal 
entries in the Hermite form of A, solving the system PA = G with the correct 
diagonal degrees gives the matrices U and H such that UA — H where H is the 
Hermite form of A. 

Theorem 4.2. Let A e F[t][V; 5] nxn with degpA^ < d and deg t Aij < e for 
1 < z, J < ri. Then we can compute the Hermite form H G F(t)[T>; S] of A, and a 
unimodular U G F[t] [£>;<!>] such that UA = H, with O((n 10 d 3 + n 7 d 2 e) log(nd)) 
operations in F 



Proof. Lemma 14. II and the following discussion, above shows that computing U 
and H is reduced to solving 0(n \og(nd)) systems of linear equations over F(t), 
each of which is m x m for m — 0(n 3 d) and in which the entries have degree 
e. Using standard linear algebra this can be solv ed with 0(m 4 e) operations in 
F, since any solution has degree at most me (see von zur Gathen and GerhardI 
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( 2003 )). A somewhat better strategy is to use the i-adic lifting approach of DixonI 



1982f ). which would require 0(m 3 + m 2 e) operations in F for each system, giving 
a total cost of O((n 10 d 3 + n 7 d 2 e) log(nd)) operations in F. □ 

As noted above, it is expected that we can bring this cost down through 
a smaller system similar to that of IStoriohannl (| 19941 ) . to a cost of 0((n 7 d 2 + 



n 5 d 2 e) log(nd)). Nonetheless, the algorithm as it is stated achieves a guaranteed 
polynomial-time solution. 

It is often the case that we are considering differential systems over Q(t)[D; S], 
where we must contend with growth in coefficients in T>, t and in the size of 
the rational coefficients. However, once again we may employ the fact that the 
Hermite form and unimodular transformation matrix are solutions of a linear 
system over Q[f]. For convenience, we can assume in fact that our input is in 
Z[t][D; S] nxn (since the rational matrix to eliminate denominators is unimodular 
in Q(t)[D; 6]). There is some amount of extra coefficient growth when going from 
A to A; namely we take up to nd derivatives, introducing a multiplicative con- 
stant of size around min((nd)\, el). In terms of the bit-length of the coefficients, 
this incurs a multiplicative blow-up of only 0(£\og(£)) where £ = min(nci, e). It 
follows that we can find the Hermite form of A G Q(t)[D; <5]™ xrl in time polyno- 
mial in n, &eg t Aij, deg-pAjj, and log||Ay ||, the maximum coefficient length in 
an entry, for 1 < i, j < n. A modular algorithm, for example along the lines of 



Li and Nemesl (jl997l) . would improve performance considerably, as might p-adic 



solvers and a more careful construction of the linear system. 



5 Conclusions and Future Work 

We have shown that the problem of computing the Hermite form of a matrix over 
F(t)[D; 5] can be accomplished in polynomial time. Moreover, our algorithm will 
also control growth in coefficient bit-length when F = Q. We have also shown 
that the degree bounds on Hermite forms in the differential ring are very similar 
to the regular polynomial case. From a practical point of view our method is 
still expensive. Our next work will be to investigate more efficient algorithms. 
We have suggested ways to compress the system of equations and to employ 
structured matrix techniques. Also, the use of randomization has been shown to 
be highly beneficial over F[t], and should be investigated in this domain. Finally, 
our approach should be applicable to difference polynomials and more general 
Ore polynomial rings. 
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